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correspond to a transition from the resonance state to the virtual state while the Coulomb forces 
gradually decrease to zero. 
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I. INTRODUCTION 



The vertex constants (VC) Gf (I is the orbital momentum, the other quantum numbers are omitted) for a •<-> b + c 
virtual decay (synthesis) and the related residues of the partial elastic scattering amplitude in the pole of a bound or 
resonance state (when a decay is real) are considered to be independent physical quantities. They play an important 
role in the theory of nuclear reactions using Feynman diagrams and in astrophysics. For the first time the vertex 
parts of diagrams, defining an input of the specific mechanism into the cross-section of a direct nuclear reaction, were 
considered in the frame of the dispersion nonrelativistic approach proposed by I. S. Shapiro more than 45 years ago 
1] . The main aim of this approach was giving up the perturbation theory, on which the distorted wave method is 
based, for processes with a strong interaction. The selection criterion for the most important diagrams is closeness 
to the observable physical region of the singularity of the corresponding amplitude along with the diagram "weight" 
produced by the product of the vertex parts. Vertex constants are directly related to the asymptotic normalization 
coefficient (ANC) of the wave function, determining (along with the "binding" energy) the system a at large distances 
between particles b and c. The Coulomb repulsion effects must be taken into account for the charged particles. In the 
present study the problem is solved in the frame of an effective-range theory. The preliminary results for the s-wave 
were published in our works 0, Q . 

As far as we know, this method of finding VC has not been considered in detail in the literature with the exception 
of [H where the authors attempt to find a relationship between the vertex constant Gf c (the upper indexes point to 
the nuclear and the Coulomb interactions) for the bound state of the nucleus a and the parameters of the standard 
effective-range approximation for the elastic scattering of charged particles b with c at arbitrary values of the orbital 
angular momentum I. Besides the binding energy £f, c of the nucleus a the parameter list includes the "scattering 
length" ai and the "effective range" r/. (We put these in quotes because the dimensions of a/ and n depend on I, 
the length dimension this parameters have only at I — 0) . Unfortunately, a serious error of fundamental importance 
was made in [4]: Eq. (25) in Q (the numbering of the formulas in this and in the next section corresponds strictly 
to the numbering in (ij), which relates the binding energy to the "scattering length" and to the "effective range", 
was written without allowance for the Coulomb interaction (!). Equation (25) which is inappropriate in the case of 
charged particles, was used there to derive a formula for the elastic-scattering amplitude {see Eq. (23) in [1]}, and the 
correct expression for cot 5f iC (E) from ||, which takes into account both the Coulomb and nuclear interactions, was 
substituted into it. As a result, the "hybrid" amplitude obtained thereby has neither a pole for the the without the 
Coulomb interaction nor for the proper binding energy value eJJ. with the Coulomb interaction taken into account. 
Nevertheless the derivative of the numerator was taken at the binding energy E cra = —the (30) in order to obtain Gf c . 
Thus the expressions (30), (31) in Q for the vertex constant of the virtual decay a — > b + c into two charged fragments 
are invalid. At first glance the situation can be remedied by using the experimental value of the binding energy 
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the in formula (31). But it is known that the parameters of the effective-range function change when the Coulomb 
interaction is taken into account. That is why one should not use formula (23) for neutral particles, which gives ef, c 
in terms of the parameters of the effective-range theory for neutral particles to find these parameters. Finding ai 
and n by using data on the vertex constant is just one of the objectives of this means using formulas (25) and 
(31) (or (34) for VC). Invalid formulas (25) and (34) were applied in paragraph 4.2 in [4j for the estimation of the 
scattering lengths and corresponding effective radii for at scattering in the p-wave by fitting the experimental ANC 
values for the channel spin Sb =3/2 and 1/2. Because of the error mentioned above the results of this paragraph are 
questionable. 

We derive formulas for the renormalized scattering amplitude and for the corresponding vertex constant in terms 
of the effective-range approximation with the Coulomb interaction taken into account. A more detailed analysis of 
the mistakes in 0] is given below. 

The prime advantage of the approach proposed by us is the possibility of finding Gf using only experimental data, 
which allows us to obtain parameters of the effective-range function Ki(k 2 ). For the arbitrary I values the standard 
effective-range expansion has the following form (see, for example, [f|) 

Ki(k 2 ) = —1/ai + nfc 2 /2 + • • • , (1) 

where the dimensions are as follows: [a;] = L 2l+1 , [rj] = L~ 2l+1 (L - is the length). The following notations are 
adopted: k = {2^,E c - m ) 1 ^ 2 is the wave number, /i is the reduced mass, E cm is the energy in the center-of-mass (cm) 
frame. For brevity below we omit the indexes of the energy in the cm frame E = E cm . The energy in the laboratory 
frame we denote as -Eq a b- Here and below we use the system of units where h=c=l. 
For the Nd system the effective-range function has the form 

K (k 2 ) = (-l/oo + C 2 k 2 + C 4 fc 4 )/(1 + k 2 /K 2 ). (2) 

The pole of the effective-range function Ko(k 2 ) in ([2|) at k 2 = —kq 2 , i-e. at the energy E cm = —Eq, where Eq = 
(3/4to)ko 2 (for the Nd system \x = (2/3)m, m is the nucleon mass), is a feature of the doublet Nd system 6-9]. For 
the nd system, the pole is at a negative energy near the threshold, i.e. Eq > 0. For the pd system the situation is much 
more ambiguous. Due to the Coulomb repulsion, the measurements at very low energies are complex and unreliable. 
As a result the information on the scattering length and position of the pole of K (k 2 ) obtained in different works is 
ambiguous and contradictory. In a few analyses of the phase shift it was found that Eq < 0. However, an analysis of 
the latest calculations within the three-body problem involving three-particle interaction gives Eq > 0, with a pole 
located very close to the elastic scattering threshold. 

In the absence of the Coulomb interaction, the effective-range function Ki{k 2 ) is related to the phase shift 5i{E) 
by the expression 

Ki(k 2 ) = k 2l+1 cot 6i(E). (3) 

In the presence of the Coulomb repulsion, the right-hand side of the formula is transformed in the well-known way 
(see, for example, [T3|) [42| : 



K^k 2 ) = k 2l+1 {c h )- 1 [[2^ 7 /(exp(2^ 7 ) - l)](coaf (£) - i) + 2 7 if( 7 )] , (4) 

where 

ff( 7 ) = *(i 7 ) + (2i 7 )- 1 -ln(i 7 ), (5) 
i 

(^r^n^+^V), m~i = i, (6) 



5^(E) = 5 l (E)-a l (E). (7) 

The following designations are used here: Si(E) is the phase shift for the sum of the Coulomb and nuclear potentials, 
<Ji(E) is the Coulomb phase shift, determined by the relation 



exp{2iai) = T(Z + 1 + i~f)/T(l + 1 - i 7 ), 



(8) 
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^fry) is the psi- function (logarithmic derivative of the gamma- function) , 7 = X/k is the Sommerfeld Coulomb pa- 
rameter, A = fiaZfrZc, a = e 2 /he is the fine-structure constant; and Zj is the charge number of the nucleus /. We 
also use the Bohr radius 0,3 — 1/A. This notation was accepted in [TTJ] . Using the explicit expression for ^(17) in the 
form of infinite sum, we can write Zf (7) as (see, for example, 

00 1 

H (l) = 77T L T-^ +1 2 Y -T-5 - ln (7) - C, (9) 

w/ exp(2?r7) - 1 ' ^n(^ 2 +7 2 ) 

where C ss 0.5772 is the Euler constant. 

Kok [ll[ was the first to show that in the presence of the Coulomb repulsion the scattering amplitude pole in the 
s-wave for the virtual (antibound) state shifts from the negative imaginary axis (link < 0) to the fourth quadrant 
of the complex plane and becomes a resonance pole at k = k Tes = Rcfc ros + ilmfc rcs . Since H (7) in ([9]) is an explicit 
function of the argument ik, another pole arises at k = — Rcfc rcs + ilmfc rcs in the third quadrant of the complex plane, 
which is symmetric relative to the imaginary axis with respect to fc res . For the singlet Np and the doublet Nd systems 
both poles are below the threshold (Re.E ros < 0), because | Imfc ros /Refc ros |> 1. The mirror symmetry of zeros and 
poles with respect to the imaginary momentum axis stems from the general symmetry properties of the S-matrix (see, 
for example, HQ)- 

The expression for cot 5q (E) , valid for real positive energies E in the physical region, is given in the textbook by 
Landau and Lifshitz (Eq.(136.11) [HI), where an expansion of type (1) is used. 

For the sake of completeness, we present the corresponding formula taken from which is valid at arbitrary I 
value and can be used for the analytical continuation into the complex momentum plane. 

Ki{k 2 ) = fc^+^c^)" 1 [[27T7/[exp(27T7) - 1]] cot Sf (E) + 7 [*(1 + i-y) + *(1 - i 7 ) - 2 In 7]] . (10) 

Eq. (fTU]) is suitable for a parametrization of the scattering phase shift. The corresponding expression for the p-wave 
is given in the work [l6| , where it was used for finding the S-matrix poles for Not scattering near the elastic threshold. 
We recalculate these pole positions more accurately in the present work using the results of the phase shift analyses 
of the experimental data known from the literature. After that we calculate the renormalized VC, which were not 
considered in [16fl. They can be used for finding the ANC of the corresponding Gamow wave functions. 

In our work [2| we study the resonant subthreshold state of the nucleus 2 He and the resonance for the aa system 



which has an extremely narrow width. The resonance positions for these states were calculated by Kok 11] using 
expansion (1). However, the vertex constants were not considered. The parameters for the pp scattering were borrowed 
from (17} . 

In works 0, 9 we also studied the 3 He bound state and the subthreshold resonances in pd scattering for the most 
reasonable sets of constants in the expansion (2) which were fitted recently in [l8l l28ll using the doublet phase shift for 
the s-wave which was calculated for few energy values in the three-body approach [2Ct l2lj . for the NN potential AV18 
and the three-particle interaction UR-IX. The necessity of using the calculation results for finding the parameters of 
Ko(k 2 ) is due to the fact that the experimental data for the pd scattering show large errors which increase while the 
energy decreases, and because the results of the phase shift analyses published in the literature are contradictory. 

In the present work the trajectories for a transition from the resonant subthreshold state to the virtual (antibound) 
state for the Np and the Nd systems were constructed with a gradual decrease of the nucleon charge (the preliminary 
calculations were fulfilled in [3] ) . As a result we demonstrate the general physical nature of those corresponding states 
which differ only in the Coulomb interaction. A comparison of the trajectories for the different systems allows us to 
understand why the role of the resonance (or the virtual level) in Nd scattering is less important than in the case of 
Np scattering. The resonance momentum trajectory at the transition from pp to np system was previously calculated 
in [22| using the Eikemeier-Hakenbroich NN potential, which is independent of nucleon charge properties, and the 
Coulomb potential multiplied by the coefficient £, which changes from to 1. The analytical continuation of the 
S-matrix onto the lower half-plane of the complex momentum was fulfilled by the Schrodinger equation solution. The 
resulting trajectory occurs close to the linear. 

II. METHOD FOR CALCULATING THE RENORMALIZED VERTEX CONSTANT 

The scattering amplitude for charged particles can be written as the sum of the pure Coulomb amplitude /coui(k) 
and the amplitude originated from a short-range (nuclear) interaction in the presence of the Coulomb field (see, for 
example, |23| ) 

/(k) = /coui(k)+/ N uci(k), (11) 
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°° 21 4- 1 

/coui(k) = J] -^r_(exp(2ia I ) - l)Pi(cos9), (12) 



2^ + 1 

/Nud(k) = ~YT exp(2^)(exp(2Mp) - l)fl(cos0). (13) 



Equation (|13[) is written for spinless particles. The generalization to the case of spins taken into account is elementary 
because the Coulomb interaction does not depend on the spin values. We consider amplitudes which are diagonal 
relative to spin and to orbital angular momentum values. A cross-section includes non-diagonal amplitudes as well 
when the spins are taken into account. In any case VC depends on one set of the quantum numbers I and s. 

The Coulomb- nuclear partial scattering amplitude with the orbital angular momentum I, which has the poles we 
take interest in, can be written as (see (fT3")0 [43j : 

/ ; c -exp(2^)/ ; , (14) 

fi = (exp(2^f - l)/2ik = l/(fccot<5f - ik). (15) 
Let us find {kcotSf — ik) from ([J]) and insert the result into (|T4"]) . (|T5"j). We obtain the following formulas: 

f?=f? N k 2l <f>i{k), (16) 

where 

fp N = me) - 2fc 2i+i ( Q7 r i 7J ff( 7 )r\ a?) 

4>i{k) = exp(2 2( r / )[2^7/[exp(2^7) - l]]^)" 1 . (18) 

Let us simplify the expression for 4>i{k). After simple transformations using the relationship (see (3.2) on page 54 
in 0) 

T(l + 1 + t 7 )r(l + 1 - t 7 ) = 7rP; + 1 ^ , (19) 

7 sn(7r 7 ) 

where P n (not the Legendre polynomial) is defined by formulas 

i 

Pi(z 7 ) = 7 2 , Pi+i(H)=l 2 U(n 2 +l 2 ), (20) 



n=l 



we receive the expression 



>iii= > rQ + i + il) y e _ wy (2i) 



The function 4>i(k), which does not contain information about the nuclear interaction, coincides with the factor at 
the renormalized scattering amplitude given in the review [25[ (see Eq.(3) in (25[). 

Correspondingly, we can write the analog of the renormalized partial amplitude for a particle scattering by a short- 
range potential in the presence of the Coulomb interaction in terms of the effective-range function Ki(k 2 ) as [in the 
following, we use a tilde sign (~) above the designation of renormalized quantities] 

k 21 

fw(k) = f? N k 21 = [Ki{k 2 ) _ 2Xk 2i H{l){cij yiy ( 22 ) 



According to this formula, it is obvious that the zero of the denominator in (|22p when k = p is an amplitude pole 
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For the standard effective-range expansion the pole position is derived from the equation (with (p} taken into 
account) 

- 1/aj + r lP 2 /2 + ■■■ = 2Xp 2l H(X/p)[{c h )- 1 ]\ k=p . (23) 

In the absence of the Coulomb interaction the pole position is defined by formula ([3]) at coti5p = i (see (15)), which 
leads to the following simple equation for the pole position in the effective-range approach with taking the first two 
terms in ([T]): 

K^p 2 ) = -l/o, + r lP 2 /2 = ip 2l+1 . (24) 



For the bound state a of two particles when one of them is chargeless this implies the equation (p — in, k — ^/2fj,eb c > 
0, €bc is the binding energy) 



l/a l +r l K 2 /2 = (-l) l K 2l +\ (25) 



which differs essentially from the equation (|23|) for charged particles. 

Nevertheless, just equation ([23)1 was used (see (25) in [J]) where in addition the factor (—1)' was lost. Correspond- 
ingly, equation (24) in [4| is also incorrect where formula (25) was used to write the scattering amplitude. In our paper 
we use the right equation (|23p from which the main equation of [11| was also derived for determining amplitude poles 
in the momentum complex plane, that is, on the nonphysical sheet of energy. The authors of [4| do not discuss the 
introduction of the renormalized VC which is real for a bound state (see (25J). In particular, the factor exp(2icr/) was 
not included in their formula (19) for the Coulomb-nuclear amplitude. We note that it is precisely the renormalized 
VC is related to the ANC by the simple relationship (se e (26|V 

Using the VC definitions which is well known from |26j | we obtain the following expression for the renormalized VC: 

G 2 = -(tt/^ 2 ) lim (fc 2 + n 2 )f lN {k). (26) 

k— Hk 

The equation (f^5|) can be rewritten as 

c 2 = (-WmV' +1 (27] 

1 ±[K l m~2\k»H( 1 )(c h )-\ = ; [ > 

where for the standard expansion ([1]) 4^ [Ki(k 2 )]k =p = r\p + • • ■ , p = in is the position of the pole amplitude for a 
bound (k > is realtor a resonance state (Imp < 0, Ren < 0). 

We considered in Q both the conventional expansion of the effective-range function K (k 2 ) in powers of k 2 (do is 
the scattering length, r$ is the effective radius, P, Q are the shape parameters), 

K (k 2 ) = -1/ao + r fc 2 /2 - Pr 3 k 4 + Qr 5 k e , (28) 

which leads to 

~ 2 2'kk/ijl 2 
= W) - [r K + 4P(r oK )3 + 6Q(r oK ) 5 + ■■■]' (29) 

and the function with a pole ^ (Co = —1/ao), when we obtain [44{ 



ip(x)~2KK 2 [~Co + C 2 K 2 -C i K 2 (2n 2 -n 2 )]/(n 2 -K 2 ) 2 ' 



The function tp(x) has the form (x = \/n = 1/obk) 

(p(x) = -1 - 2x + 2xH'(x). (31) 

Let us remember that due to (15) the positions of the poles (fc rcs = in) of the scattering amplitude /(fc) for the 
bound or resonance (or virtual) states are determined by the condition 

cat8?(E) = i. (32) 
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III. VERTEX CONSTANTS FOR PARTICULAR NUCLEI 

The case of the s-wave (/ = 0) was considered in our papers [HIHj where all the necessary formulas were given as 
well as the numerical calculation results for the nuclear systems NN, Nd and oca. Several results from [2, 3] are given 
below for the sake of completeness. In the present paper we study more thoroughly the trajectories of the resonance 
momentum and the renormalized VC while the charge of the one of the particles gradually goes to zero. 



A. The proton-proton system. The ground state of the 2 He nucleus 

The results of the transcendental equation solution Q23|l with the parameters set, taken from [l7| (without taking 
inaccuracy into account), found in 0, are fc2 Hfi = (0.0644 - iO.0871) fm -1 , i?2 Hc = (-142 - i465) keV. They differ 
only in the last figure from those obtained in [ll]: fc 2Ho = (0.0647- i0.0870) fin -1 , E2 Hc = (-140 - i467) keV. In 
for the first time the renormalized vertex constant squared was also found G 2 = — (0.060 + i0.051)fm. The calculation 
results change only weakly if the shape parameters P = Q = 0. 

Let us now consider the transition from the pp system to a system without Coulomb interaction. Either nn or np 
system can be taken from the isotopic triplet. For the sake of certainty, we choose the singlet np system, for which 
(in the approximation P = Q = 0) a% p = (-23.719 ± 0.013) fm, r£ p = (2.76 ± 0.05) fm (see []~7]). With this set of 
parameters, the virtual level energy E v = —66 keV. 

In the pp np transition was set by the linear relations 

a (£)=a^ + C(l-£); (33) 
r (C)=r^ + r7(l-C). (34) 

This approximation leads to a highly nonlinear trajectory of the resonance pole fc res (£) in the momentum complex 
plane, in contrast to the result of [23 |. It is known (see, for example, [l7j). that the effective radius only changes 
slightly when the Coulomb interaction is taken into account; therefore (l34l) is a good approximation. The decisive 
factor is the dependence of the scattering length on the product of charges, which significantly differs from the linear. 

A monograph by Brown and Jackson [17| gives an approximate formula for the relation between the pp and np 
scattering lengths derived in accordance with the calculations performed long ago by Landau and Smorodinsky (see 
(4.29) and reference in [l7|). In [l7| it is noted, that, in spite of crude approximations, this formula with an accuracy 
within a few percent is valid for the potentials, which permit us to describe well the pp scattering data. Analogous 
formula had been derived in the work by Schwinger [27| . 

Introducing the coefficient £ into this formula by the replacement A — > £A, we obtain the following relation, 
determining the dependence of the scattering length a p (£) on £ Q: 

l/a p (0 = l/a" p + 2A£ (ln(2i?A£) + 0.33). (35) 

The value of R was not found in [171 ] . It was only noted that the dependence on this parameter should be weak 
because R is under the logarithm sign. Solving the transcendental equation (|35|) at £ — 1 we find R = 1.74 fm. This 
value is within the margins of the radius spread for different potentials with two parameters given in (l7| for the 
singlet np system. The dependence in the form (|35|) is approximated well by the relation [2j 

l/a p {£,) = {l-e)/a7 + e/al p (36) 

at (3 = 0.7. Clearly, nonlinear dependence %,(£) is shown in FigfTJ 

Finding out the dependence of the parameters of the effective-range function on £ is a dynamic problem. Let us 
introduce into the Schrodinger equation an interaction £Vc(r), instead of the Coulomb potential Vc(r), where the 
coefficient £ ranges within < £ < 1. Thus, one must replace the Sommerfeld parameter A — ¥ £A in all formulas. We 
obtain the results shown in Fig. [T]by solving the Schrodinger equation in the continuum for the Np and Nd scattering 
to find the energy dependence of the effective-range functions and their parameters. We use the Yukawa potential 

V(r) = ~V (R/r) exp(-r/R), (37) 

as a nuclear interaction for which the parameters of the neutron-proton interaction in the singlet state are given in 

Let us now construct the trajectories for k TCS (£,) and G 2 (£). We find the parameters of the effective-range function by 
fitting the phase shift values calculated at low energies where the effective-range expansion is valid. The corresponding 
pole trajectories in the complex momentum plane are given in Fig. [5] They are close to the linear as in (22j . 
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The calculation results show that, in accordance with the conclusion of |17[, it is only the dependence of the 
scattering length a p (£) on £ that is important. The effective radius changes only slightly and can be described by 
the linear dependence (|34|) . The corresponding curve is situated a little lower than the curve following from the 
Landau-Smorodinsky-Schwinger (j55| (or at /3 = 0.7), because in this case the pp scattering length differs from 
the experimental one. 

The trajectories G 2 (£) in the complex plane of the renormalizcd vertex constant for the Np system given in Fig. 
|3] are also close to the linear. They are symmetric with respect to the real axis - that is, a complex conjugate. 
The contribution of the poles, located symmetrically with respect to the imaginary momentum axis, to the partial 
amplitude on the real energy axis is real in this case, as it should be for real interaction potentials. 



B. The doublet state of the pd system in the s-wave 

With the strong correlations between the parameters C2, C4, Eq ([2]) and the scattering length ao taken into account 
(see, for example, 0, UMl) it is suitable to write down the effective-range function in the form 

„ ,, 2 , (-1 + c 2 fc 2 + C4 fc 4 ) -1 + b^ + hE 2 

Kaik } = (a + d 2 ki) = (oo + iiiO ■ (38) 

The parameters which appear in this formula have a much weaker dependence on ao- An approximation of the results 
of three-body calculation of the phase shifts in the explicit form ([2]) and (|38|) can be considered as an analytical 
continuation of the partial scattering amplitude onto complex energy E values. This allows us to find the whole set 




FIG. 1: The ao(£) dependencies of the singlet Np scattering length (two bottom curves) and of the doublet Nd scattering length 
(top curve, triangles) on the coefficient £ (multiplier before of the Coulomb potential). The results of the ao(£) calculations are 
obtained using (35) (following the Landau-Smorodinsky-Schwinger relation between and ctg P ) (circles) and by solving the 
Schrodinger equation with the Yukawa potential (squares). 
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of the low energy characteristics for the Nd system. The set includes the positions of the scattering amplitude poles 
for the bound and for the resonance (or virtual) state which are defined by zeroes of the denominator of (22) due to 
the condition l|32[). along with the pole position (—Eq) of the function K(k 2 ). 

The corresponding residues at the poles are connected with the nuclear vertex constants G\ and G 2 (see [2H, [2r|) 
or with the asymptotic normalized coefficients C 2 and C 2 (see, for example, 29]) for the nuclei 3 H and 3 He, which for 
the sake of brevity are denoted by the letters t and r correspondingly. As is noted in the introduction, the available 
phase shift analyses of the experimental data for the pd scattering do not give an unambiguous result. In the first 
place, this affects the scattering length. Because of this, the analyses of both the modern three-body calculations 
with three-particle forces taken into account [2(| [U [3(| and the calculation results in the frame of N/D method [3~fl ] 
were fulfilled in QS EH- 

Different physical observables were used for the parameter fitting: the s-wave doublet phase shift for the pd 
scattering at the smallest energy values considered in the literature, the binding energy of 3 He, and the position of 
the subthreshold resonant pole, corresponding to the virtual state of 3 H, which was calculated previously by the N/D 
method in [3l|, where the residue of the scattering amplitude was calculated as well. The different fitting variants 
were given in the tables of [HI, EH . 

The renormalized VC values for the selected fitting variants are given in Table 1 of our work Q where the right 
-E rcs values and the parameters of the effective-range function were presented as well [45] . In the second block of this 
table the variants for the Argonne NN potential (AV18) with three-nucleon forces taken into account in the Urbana 
form (UR-IX) (see References in [2(], [2l|) are considered for the value 2 a p d = 0.024 fm, obtained in [3(| published 
later where the pole position of the effective-range function is also found (Eq — 3.13 MeV). This simplifies the fitting 
of the remaining parameters. In this case it is possible to reproduce precisely the scattering length value. Meanwhile 
the approximate phase shift behavior at low energy well describes the original three-body results for the interaction 
AV18+(UR-IX) no matter whether the binding energy is fitted or not. For this reason, we consider 2 a p d — 0.024 fm 
as the most reliable theoretical estimation to date. 

The resulting parameter value C2 = 2 a p dC2 — —58.1 fm 2 found by us can be compared with the result of the analysis 
of the phase shift calculated in the energy region E cm < 450 keV by the Pisa group (see [30]) using the formula (|3"5)) 
when C4 = 0. The result obtained which corresponds to the value C2=— 56.7 fm 2 agrees well with the results of the 
analysis in [ijj], especially for the variant No 10 (c2 = —56.6 fm 2 ) (see Table 1 in [2]). 




FIG. 2: Trajectories of fc rcs for the transition from the resonance pole to the virtual (antibound) state at a gradual decrease 
in the Coulomb interaction to zero. The top curves are trajectories for the Np system. The points correspond to different £ 
values in the interval 0—1. The circles present the calculations, obtained by formula (|35[) ; the squares for the Np and points 
(bottom curves) for the Nd systems are obtained by solving the Schrodinger equation with the Yukawa potential. The thin 
straight lines are the angle bisectors of the third and fourth quadrants of the complex plane. 
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1. The bound state 3 He characteristics 

The renormalized VC for the 3 He bound state calculated in in the effective-range approximation (2) occurs 
to be greatly underestimated (by about a factor of 2) in comparison with the results of both experimental analyses 
and theoretical calculations. This is due to the fact that the 3 He bound state pole is situated far away from the 
convergence region for the expansion of the numerator ([2]) in the powers of k 2 . 

Unfortunately, the majority of available publications on the VC for 3 H and 3 He are some 20 years old. In addition, 
different papers report values of either the renormalized VC G 2 or the asymptotic normalization constant (ANC) 
C 2 . Table 1 in 2] presents both the renormalized VC and the ANC for the 3 He bound state. The following relation 
between the VC and the ANC (see, for example, [28, 29]) is used: 



G 2 = (27/A)ttk t X n 2 C t 2 , (39) 

where XN—h/mc, k t is the wave number of the bound 3 He nucleus. This relation corresponds to the definition of the 
asymptotic behavior of the normalized-to-unity relative motion wave function of the proton and the deuteron in the 
form (for the s-wave) 

V(r) -»■ C T N ZR [W_ lil/2 {2K T r)] /r, r -> oo, N ZR = (2k t ) 1/2 , (40) 

where r is the relative distance (Jacobi coordinate) between the nucleon and the center of mass of the deuteron. 
This definition of the ANC corresponds to the value C T = l in the zero range approximation, when the wave function 
coincides with its asymptotic form (|40[) for any ranges down to r = 0. 

For comparison, we give the VC and ANC data from [2{|, where the ANC is introduced by the relation 

il>(r) -»■ C T N w [W^ 1/2 (2Krr)}/r, r ^ oo, (41) 



(S Q 




-0.06 - 



-0.06 -0.04 -0.02 0.00 

Re(Gi) (fm) 



FIG. 3: Trajectories as in Fig. [2] for the square of the renormalized vertex constant, G 2 en = (G) 2 , calculated with the Yukawa 
potential for the Np (squares) and Nd systems (on the left and right, respectively); the trajectory 5(G) 2 is also shown for the 
Nd system. 
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TABLE I: The nucleon separation energies B t and B T (in MeV) for the nuclei 3 H and 3 He, the vertex constants G\ and 
renormalized G T (in fm), the asymptotic normalized coefficients squared Cf and calculated in the two-body model for the 
different effective nd potentials (for the Hulthen potential the values are more accurate than in comparison with those given in 
0)- 



Potential 


B t 


c? 


Gt 


B T 




m 
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TJTn 1 + Vi nil 

nuitnen 
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Q Al 
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1 Ad 
1.40 


O.O / 


o.Do 


1 A « 
1.40 


U.00 


Yukawa 


6.26 


3.06 


1.28 


5.64 


3.28 


1.31 


0.66 


Tomio et al. (B) 


6.26 


3.56 


1.49 


5.60 


3.77 


1.50 


0.65 


Tomio et al. (A) 


6.26 


3.94 


1.65 


5.58 


4.18 


1.66 


0.68 


where |46| 




(N w ) 


OO 

--JV 


- 1 ,i/2{^, T r)] 2 dr. 






(42) 







The analytical expression for Nw is given in HH, where N w /N Z r = 1.024 for 3 He (see This factor is taken 

into account below in the calculation of the VC from the results for C 2 from (|4"Tj) . The calculations in [29] were 
carried out by solving Faddeev type equations in the configuration space for the Malfliet-Tjon (I III) (MT) and Reid 
soft-core (RSC) NN potentials. The results are as follows (e T is the binding energy of the 3 He nucleus in MeV, G 2 
in fm): 

(MT) e T = 7.87; C 2 = 3.90; G 2 T = 1.63; (43) 

(RSC) e T = 6.39; C 2 = 3.14; G 2 = 1.12. (44) 

Coulomb effects are calculated in (|43| and (|44|) within the point charge approximation. Allowance for spread of the 
charge does not affect the VC and ANC values to a precision of one-hundredth. The results for the RSC potential 
are considerably underestimated in comparison with those for the MT potential, first of all because the nucleus 
is underbound for the RSC potential. In [30j the binding energies of the 3 H and 3 He nuclei which are near the 
experimental values were found for the NN potential AV18 with the three-particle forces UR-IX taken into account. 
We give below the results (the units are the same as in (1431) and flSJ}) for 3 H and 3 He with the interaction AV18+ (UR- 
IX): 

e t = 8.49; C 2 = 3.44; G 2 t = 1.44; e T = 7.75; C 2 = 3.53; G 2 . = 1.39. (45) 

It has long been established in the literature that there is a correlation between the binding energy and the vertex 
constant (see, for example, [HI and references therein). Therefore, a good reproduction of the experimental e T value 
is important for obtaining a reliable value of the VC. The two-body potential model [32|,[33j], where the binding energy 
and the scattering length are adjustable parameters, meets this condition. 

Recently we have carried out the corresponding calculations for various potentials (3~H , where the potential param- 
eters were given. The parameters of the Hulthen and Yukawa effective nd potentials were fitted to the experimental 
values of the 3 H binding energy and to the doublet nd scattering length 2 a nd . On the assumption of charge indepen- 
dence, the same nuclear potential was used for the doublet pd system. In addition, calculations were carried out for 
the potentials proposed in paper by Tomio et al. [33| (versions A and B). In the case of version B, we refined the 
parameters beforehand to get the correct doublet nd scattering length (see reference in [34| ). 

The results of the calculations are presented in Table U [47[ . It is noteworthy that G\ and G 2 are close in value 
despite the noticeable difference of the ANCs. The Coulomb difference of the 3 He and 3 H binding energies (« 0.7 
MeV) is in rather good agreement with the experimental value A£' c =0.76 MeV. 

The results for 3 He agree with the conclusion in [TH, [l9], 35] concerning the convergence domain of expansion (2) 
limited by the energies |-E C m| < 0.74 MeV. This domain is determined by the position of the nearby singularity of 
the partial wave amplitude for the pd scattering at the energy E cm — —6^/3 (e^ is the deuteron binding energy) 
corresponding to the Feynman diagram for the one-nucleon exchange. The poles for the bound states of the Nd 
system turn out to be far beyond the convergence domain. At the same time the resonance of 3 He lies within the 
convergence domain as the pole for the virtual triton and therefore its position and characteristics should be less 
sensitive to the variant of fitting and thus more reliable. We think that the best are the fits No 7, 8 in p~9j ] . 



11 



2. Subthreshold resonance of the pd system 



Our calculation of the energy of the subthreshold resonance for the the parameter set of the effective-range function 
taken from [3(| leads to the value Ef^ s = —(0.315 + zO. 102) MeV, which differs only in the third decimal figure from 
the corresponding value Eg* = -(0.319 + i0.099) MeV for the variant No 8 in Table 1 of 0, where the binding energy 
of 3 He serves as an additional fit value. Let us compare our result Ef^ —(0.32 + «0. 10) MeV with other results 
published in the literature. 

In [3l| the value E$* = -(0.432 + i0.032) MeV was obtained by the N/D method. The absolute value of the real 
part is slightly overestimated in the N/D method whereas the absolute value of the imaginary part is smaller by 
about a factor of three. In H| the value Ef^ s = -(0.432 + i0.56) MeV is given, which was obtained by solving the 
Faddeev equation for the Eikemeier-Hackenbroich NN potential. Its real part coincides with the result of [3l|, but 
the absolute value of the imagina ry p art is larger than in [19J by about a factor of five and is an order of magnitude 
larger than in the N/D method [31j, which calls for explanation. Note that the position of the virtual pole of the 
triton found in [22| (B v = 1.62 MeV) is also considerably larger (by a factor of more than 3) than other estimates 
reported in the literature. 

Let us now construct the trajectories of the position of the pole and residue in it at the transition from pd to nd. 
As in the case of the nucleon-nucleon system, we use the Yukawa-potential model (j3"7| with the parameters and the 
correponding values of the physical quantities ANC and VC calculated by solving the Schrodinger equation for the 
bound and continuum states (lengths in fm, energies in MeV): 

(Yu,nd) R = 2.77, V = 17.46; C t 2 = 3.06, G? = 1.28. (46) 

The parameters in (|4"6")l are fitted to the triton binding energy (e* = 8.48 MeV) and to the doublet scattering length 
( 2 Ond = 0.65 fm) for the nd scattering in the s-wave. This model describes well the energy behavior of the phase shift 



for the doublet Nd scattering [34| and other low-energy characteristics which were discussed above. 

We introduce the coefficient £ into the product of charges and find the dependence of the s-wave phase shift 6q(E) 
on E in the range of the convergence of the effective-range function Kq (k 2 ) at different £ values (0 < £ < 1). 

In Fig. [4] we depict the dependencies K^^Ei^) on the energy in the range < 0.3 MeV, where the Coulomb 
interaction effects are most pronounced. The calculations were performed with the set £ = 0, 0.2, 0.4, 0.6, 0.8 and 
1.0; these values correspond to the downward sequence of thin curves in Fig. 4 (the top and bottom curves are for 
nd and pd scattering, respectively). The points on the curve correspond to the phase shift values found by solving 
the Schrodinger equation in a continuum with the Coulomb potential £Vc(r). The curves connecting the points are 
the results of fitting using formula We took the set of fitting parameters c 2 = aC 2 ,c 4 = aC*4 and E /a, which 

depends weakly on the scattering length a. It was found that the fitted aC^ values are unstable: they change abruptly 
in magnitude with a variation in £ taking large values due to the smallness of fc 4 in the range of E under consideration, 
and even change their sign. This behavior is in agreement with the conclusion of [l|| according to which the parameter 
C4 cannot reliably be found without fitting the binding energy of the Nd system. In the range Ei a \, < 0.3 MeV, one 
can assume that aC±=0. Such an assumption was made in [30l ] in determining the Kq(E) parameters for the pd 

system. To fit the other parameters a, aC'2 and Eo/a it is sufficient to use Kq (E) at low energies in the region 



-Blab < 0.1 MeV where the convergence of the expansion is secured. It can be seen that the found parameters for 
the considered set of £(£?) make it possible to reproduce nicely the Kq^(E) values calculated from the Schrodinger 
equation at higher energies in the range under study. 

It is remarkable that the observed effect of consideration of the Coulomb interaction within the two-body model 
is in qualitative agreement with the results obtained within the three-body problem for the nd and pd systems. 
The corresponding dash-line curves are shown in Fig. 2] together with our results for Kq(E). There is very good 
agreement for the nd scattering (the top curve). The significant observable differences for the pd scattering are related 
to the difference of the close-to-zero value 2 a p d — 0.024 fm in [30| from 0.136 fm, the value obtained in this study 
within the two-body Yukawa-potential model. 

The first estimate of the doublet pd scattering length within the two-body model was made in (33[; the result 
( 2 a p d = 0.15 fm) is close to that obtained by us using the Yukawa potential. Let us remember that the potentials 
considered in _33] include three parameters fitted from the doublet nd scattering length ( 2 a n d = 0.65 fm) and the 3 H 
and 3 He binding energies. 

Although there is no reliable experimental result for 2 a p d at present, we think that the three-body theoretical result 
in [30] is currently the best estimate. 

Therefore, for the pd system one can fit the parameters of the Yukawa and Hulthen potentials by analogy with the 
fitting for the nd system: from the experimental binding energy (e T = 7.73 MeV) and the theoretical pd scattering 
length ( 2 a pd = 0.024 fm). 
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As a result of such fitting for the Yukawa potential we found the following parameter values and the corresponding 
physical quantities calculated by solving the Schrodinger equation for the bound state and continuum (lengths in fm, 
energies in MeV) 

(Yu,pd) i? = 2.82, V = 16.93; 2 a pd = 0.026, e T = 7.73, C 2 = 3.30, G 2 = 1.30. (47) 

Obviously one can ignore the difference in 2 a p d (0.024 and 0.026). The values of C T 2 and especially of the G 2 in (|47|) 
differ only a little from those reported in (l4l)|) and in Table 1 for the nd system. We also calculated the features of the 
subthreshold pd resonance for the parameter set (|47j) using the effective-range function values from the Schrodinger 
equation by this function fitting for the set of the energy E'lab in the range (0.02-0.16) MeV using the formula (|38|l . 
With the parameters obtained in this way and when = (&i = 1.216 MeV -1 , di = 5.112 fm/MeV) we found 
Eg* = = -(0.310 + i0.074) MeV and G 2 = (-0.00463 ± i0.00215) fm which is in reasonable agreement with of our 
analysis of the three-body calculations [30] made by the Pisa group of physicists: Ef^ s — —(0.319 + i0.099) MeV, 
G 2 es = (-0.00596 ±i0.00358) fm. 

The differences between the corresponding parameter pairs which resulted from the breakdown of the charge inde- 
pendence are also small, although they lead to a rather high change of the scattering length, which is known as the 
value most sensitive to parameter change. 




FIG. 4: The dependence of the doublet effective-range function Kg (E) on the energy _Ei a b for the Nd scattering. The sequence 
of the thin curves from top to bottom corresponds to the set of values £ = 0, 0.2, 0.4, 0.6, 0.8 and 1.0. These calculations are 
fulfilled for the Yukawa-potential two-body model with the parameters (|46[) of the nd system usin g th e a p proximation (f2| when 
d = 0. The thick dashed curves are obtained from from an analysis of the three-body results [2(], [21[ for the Argonne AW 
potential (AV18) taking the three-particle Urbana forces (UR-IX) into account. The top dashed curve is for the nd scattering, 
the bottom dashed curve is for the pd scattering. The lowest solid thick curve shows the calculations for the Yukawa potential 
with the parameters Q47I) for the pd system. 
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We calculate the function Kq '(E) with the parameters (l47l) at £ = 1 (see Fig. [4|. A comparison with the three- 
body curve in the same figure shows that in the energy region under study the two-body model with the Yukawa 
potential when the parameters are fitted from the three-body binding energy value and the doublet pd scattering 
length given in (|47)l leads to a larger difference between the values of the function Kq(E) at £ = and £ = 1. This 
means that the role of the Coulomb repulsion is somewhat overestimated. But the observable difference between the 
two-body and three-body curves at £ = 1 is well known to be smaller than the experimental error. In the low-energy 
region considered reliable measurements of the pd scattering actually does not exist due to the high Coulomb barrier. 

With the parameters of the effective-range function K^p (E) found in such a way for the Yukawa potential for the 
nd system we calculate fc res (£) and G 2 (£). 

The charge dependence of the scattering length 2 apfd(£.) ( the most important for the trajectory constructing) which 
we calculate within the two-body Yukawa-potential model with the parameters (|46j) is shown in Fig. [TJ This function 
2 aNd(£,) is well approximated by a quadratic polynomial 



The changing of the scattering length with the charge variation is quite sudden: the ratio of the scattering lengths 
at £ = and £ = 1 is close to 4, whereas the corresponding ratio in the case of the Np scattering is near 3. The 
changing of the reciprocal quantity [ 2 ajV£;(£)] _1 is quicker. It can be described by the exponentially increasing function 
at £ — > 1 or as a function with a pole at £ > 1. It should be stressed here that the relation (l3"5"j) is by no means 
universal. Obviously, this is not valid for the system in which the effective-range function has a pole close to the 
scattering channel threshold. However, because of the scale difference, the changing of the Nd scattering length with 
the charge variation looks much weaker in Fig. [T] compared with the dependence on £ of the Np scattering length. 
The corresponding curves differ not only quantitatively but also qualitatively. 

The corresponding trajectories of the subthreshold resonance momentum and of the renormalized vertex constant 
for a transition from the spin-doublet proton-deuteron system to the neutron-deuteron system are shown in Figs. 
[5] and |3J Explicit qualitative and quantitative differences are seen between the corresponding trajectories for the 
nucleon-nucleon and nucleon-deuteron systems. In the latter case, at a transition from the resonance to the virtual 
state, the pole trajectory is shifted not upward but downward and is more remote from the physical energy range. 
The | G 2 | magnitude for the Nd system is about factor 5 smaller than for the Np system (for clearness in Fig. 
131 the trajectory (5 G 2 ) is also shown for Nd system). The charge dependence of G 2 also differs from the case of 
nucleon-nucleon system. 

Due to the correlation between Eq and 2 a p d, the three-body Eq value also significantly differs from the two-body 
model result. There is a large difference between the values of parameter Eo/ 2 a p d which is less sensitive to a change 
in the scattering length 2 a p d- The Eo/ 2 a p d — 0.13 MeV/fm was obtained in the three-body approach [3(| whereas the 
value of Eo/ 2 a p d ~ 0.2 MeV/fm was found in the present study on the basis of the two-body model with the Yukawa 
potential. 

IV. THE p-WAVE RESONANCES IN THE SCATTERING OF THE NUCLEON BY 4 He 

The resonances in the Na scattering near the threshold have been studied in many experimental and theoretical 
works (see, for example, [36] and references therein). In [36] the results of the R-matrix analyses were given not only 
for the elastic scattering but also for the stripping and pick-up reactions. The corresponding results for the resonant 
energies and widths are also presented together with the experimental errors. It should be noted that the R-matrix 
approach used contains a considerable number of parameters. Including data on nuclear reactions in the analysis 
leads to additional errors due to its dependence on the model used to describe the reaction considered. 

As we note in the introduction, it seems that the effective-range expansion was used for the first time in [r| to find 
the pole positions on the unphysical energy sheet which correspond to the resonances in the p-wave amplitude of the 
Na elastic scattering. We limit ourselves in the present work to considering these p-wave resonances when formula 
(HJ) has the following form: 



2 a Nd (0 = 0.667 - 0.273£ - 0.244£ 2 . 



(48) 



K^k 2 ) = k 3 (l + 7 2 ) [[27T7/[exp(27T7) - 1]] cot 6 f(E) + 7 [#(1 + ry) + *(1 - i 7 ) - 2 In 7]] . 



(49) 



For the renormalized VC we correspondingly derive the expression: 




(-27T/MV 



(50) 
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TABLE II: The complex values of the resonance CM energies E r — 


iY/2 (in MeV) for f = 3/2" and f 


= 1/2" 


relative to the 


threshold in the channel a + N ( He, Li ) and of the T matrix residues Ri = 


\Ri \exp(itpi ), ((/pi in deg 


•)• 




Method Nucleus f 


E r 


L 




tpi 


Eff. rad. [16] 


0.778 


0.639 


- 


- 


N ID [371, [251 


0.697 


0.542 


0.160 


132 


(d,-y) [381 5 He 3/2" 


0.80 


0.65 


- 


- 


i?-function [39] 


0.771 


0.644 


- 


- 


i?-matrix [361 


0.80 


0.65 


- 


- 


Eff. rad. (present paper); 


0.778 


0.641 


0.171 


-133 


Eff. rad. [16] 


1.999 


4.534 


- 


- 


7V/D [37], [25] 


1.875 


5.296 


0.237 


177 


ii-function [39] 5 He 1/2" 


1.970 


0.218 


- 


- 


i?-matrix [361 


2.07 


5.57 


- 


- 


Eff. rad. (present paper) 


2.000 


4.533 


0.199 


178 


Eff. rad. [16] 


1.637 


1.292 


- 


- 


iV/£> [37], [25] 


1.655 


1.278 


0.304 


-136 


(d,7) [38J 5 Li 3/2" 


1.72 


1.28 


- 


- 


i?-matrix [361 


1.69 


1.23 






Eff. rad. (present paper); 


1.630 


1.437 


0.2/8 


1 A O 


Eff. rad. [16] 


2.858 


6.082 






jV/Z? [37], [25] 5 Li 1/2" 


2.691 


d A A Q 


0.313 


-178 


i?-matrix [361 


3.18 


6.60 






Eff. rad. (present paper); 


2.34 


p. m 


0.287 


4.2 


where the notations are the same as in [T3|: 










Q(fc) = 7 (l+7 2 )fc 3 [*( 27 ) 


+ *(-fy) 


-21n 7 ], 




(51) 


C 2 = 27T 7 (1 + 7 2 )/[exp(2vr7) - 


!]• 




(52) 



The results of our calculations are compared in Table [TT] with those obtained in [l6[ and [37| (see also (25[) as well 
as in the more recent work [38] where the experimental data on the radiation capture of the deuteron by 3 H and 3 He 
was analyzed. As far as we know, the renormalized VCs for the resonances discussed here were previously calculated 
only in [33] using the N/D method and also presented in 25 1. In Table [TT] the quantity R\ for the residue of the 
T- matrix at the pole is given as well. It is connected with the renormalized G 2 by the relationship (see (25| ) 



Gf 



(53) 



A comparison of our results with those presented in [25| for the N/ D method shows some differences in the resonance 
positions and their residues. The absolute values of the residues more or less agree with each other, while the angles 
ip almost coincide in two cases but differ in the other two. The tp values for the states 5 He (j = 1/2) and 5 Li(j = 3/2) 
are in good agreement whereas for the state 5 He (j = 3/2) the angles <p differ mainly by their sign. This means that 
the residues are a complex conjugate. But for the state 5 Li (j = 1/2) the difference between them amounts to about 
180 deg. This means that the residues differ approximately in their signs. 

Taking the Coulomb interaction into account leads to the sign change of the angles given in [25| . In our calculations, 
the corresponding change of tp more correlates with the value of the total angular momentum j whereas the magnitude 
of the residue does not depend strongly on j. The inclusion of the Coulomb interaction leads to a marked increase in 
the absolute value of the residue. Apparently, the noted differences between the tp values are not fundamental. 
It is known that resonance poles exist in pairs which are mirror symmetrical relative to the imaginary axis of the 
momentum, the corresponding residues being a complex conjugate. The sign difference of the residues also does not 
chan ge t he input of the poles into the cross section if there is no essential interference with a nonresonant amplitude. 

In [l6jj the phase shift analysis of the Not scattering presented in [4(3] was used. We compare two variants of the 
phase shift analysis of the Not scattering which differ in their maximal energy values: (A) E cm < 3 MeV for the not 



15 



scattering and E cm < 5 MeV for the pa scattering [|l| and (B) E iah < 21 MeV [ioj]. We include in Table HQ only 
our results for variant B. For the states when j = 1/2 the area of the fast decrease of the phase shift is situated at a 
higher energy (see Fig. [5]). Furthermore, the resonance width is almost ten times greater than that for the resonance 
when j — 3/2. At the same time, the phase shift does not pass the value ir/2. 

The short energy interval in variant A does not cover the area where the phase shift decreases quickly due to the 
resonance pole position. Due to this, variant A docs not permit a correct reproduction of the parameters for the 
effective-range function. When j = 3/2 the short interval is representative enough to give reasonable parameters for 
K(k 2 ) and for the corresponding resonances. The increase of the energy interval for the fitting parameters up to 
energy 21 MeV leads to slight differences in the results when j = 3/2. This testifies to the fact that the resonance 
poles considered are situated in the area where the effective-range expansion converges. 

A comparison of our results with those in [l6j shows excellent agreement for the na scattering but considerable 
disagreement for the pa scattering. We believe that our calculations are more accurate due to the progress in 
calculation methods for equation solving, including the function tp. In Table |TT] the complex resonance energies 
calculated using the R-matrix method are generally in reasonable agreement with those obtained by other methods 
(see also [36J). 

It should be noted that for the pa scattering the agreement of our results using the analytical continuation of the 
effective-range function with other results published in the literature is not as good as that between the latter and 
the results in [jfjj . 



V. CONCLUSIONS 



Thus, in this study the correct expressions are derived for the renormalized vertex constant for the decay of a 
nucleus into two charged particles in a state with the arbitrary angular orbital momentum within the effective-range 
theory. However, it should be noted that, since this theory is valid only for low enough energies when an expansion 
of the effective-range function in powers of k 2 converges, there is also a limitation on the upper value of the orbital 
angular momentum. It is likely we can largely limit ourselves to considering s and p-waves. For the s-wave ( see 
also [2] ) we consider both the standard effective-range expansion and the expansion with a pole of the effective-range 
function, which is necessary in the case of the Nd system. The formulas derived are applied to the bound state of 
3 He nucleus and to the subthreshold resonance in the pd scattering which corresponds to the virtual level of 3 H and 
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FIG. 5: The results of the phase shift analysis [4lJ of the nucleoli scattering by the nucleus 4 He for the p-wave: the solid circles 
and squares are for na scattering when j = 3/2 and j = 1/2 respectively, while the inverted and right triangles stand for pa 
scattering when j = 3/2 and j = 1/2, respectively. 
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also to the resonance of the pp system. Trajectories are constructed in the complex planes of the momentum and 
of the renormalized vertex constant (real for neutral particles) at the transitions from the subthreshold resonance 
in the ground singlet s state of the pp system to the antibound np system and from the excited doublet state of 
the pd system to the virtual triton state at a gradual decrease in the Coulomb interaction to zero. The constructed 
trajectories demonstrate the general physical nature of the corresponding states, which differ only in the Coulomb 
interaction. 

The significant differences between the corresponding trajectories for the Np and Nd system are related to both 
the larger remoteness of the virtual triton pole and to the different behavior of the scattering lengths <z(£) that depend 
more strongly on charge: with £ increasing a(£)(< 0) quickly increases for Np scattering whereas aj\rd(£)(> 0) even 
more quickly decreases for Nd scattering. Because of the different scales of variation, the change for the Nd system 
looks slow and almost a linear decrease in Fig. [TJ It is remarkable that the relatively sharp decrease in ajvd(0 with 
£ increasing due only to the Coulomb interaction, which leads to a reasonable value close to zero, was obtained by 
us within the simple two-body model with two parameters, fitted from the triton binding energy and doublet nd 
scattering length. The doublet pd scattering length found is in good agreement with the results of modern three-body 
calculations. The first result along these lines was obtained in [33] where the three-parameter two-body model was 
proposed (the third parameter was fitted from the binding energy of 3 He). The success of the two-body model in a 
qualitative reproduction of the doublet scattering length means that the internal deuteron structure is apparently of 
no importance at the nucleon and deuteron zero relative energy. 

The parameter fitting for the Yukawa potential from the three-body pd scattering length 2 a p d considerably improves 
the agreement between the two-body and three-body effective-range functions at very low energies -Ei a b < 0.2 MeV. 

The reason for the difference in the trajectories for the Np and Nd systems is that a singlet deuteron is in the 
ground state with zero spin, for which the np interaction differs from that in the spin-triplet state, leading to the bound 
deuteron. The ground state of an Np system can be described within the standard effective-range approximation. 
The virtual triton corresponds to the excited state in which all quantum numbers, except for the principal quantum 
number, are the same for the real and virtual tritons. So for the virtual state one has the same effective nuclear 
potential as for the bound state, and the function K^\k 2 ) has a pole near the threshold. 

The | G 2 | value for Nd system is smaller than for Np by a factor of about 5. Along with the larger remoteness of 
the pole position from the physical region, this fact explains why the Nd subthreshold poles are much more difficult to 
find in experiments in comparison with the well-studied Np scattering poles. The trajectory of fc r cs(£) which follows 
from the approximate formula by Landau-Smorodinsky |l7| and Schwinger (27| . relating pp and np scattering lengths, 
is in good agreement with our result for the singlet Np system in the model with the Yukawa potential. The difference 
is due to using the nuclear potential for the np system which leads to a small change of the Np scattering length at 
£ = 1 in comparison with the experimental pp scattering length. 

Finally, we recalculated the resonance positions E vcs — iT/2 in the nucleon scattering on 4 He in the p-wave within 
the effective-range approximation and found the corresponding residues at the poles which were then compared with 
the results [37,25] obtained in the N/D method. 

For the pa scattering when j = 1/2 the E res value decreased slightly and when j = 3/2 the T value increased a 
little in comparison with the results of the work [l6| where the effective-range approximation was also used. 

The expressions obtained and the results of the numeral calculations of the vertex constants for nucleus decay into 
two charged fragments can be applied to the theory of reactions using Feynman diagrams for a description of the 
process mechanism, and in the analysis of astrophysical nuclear synthesis reactions. 
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A misprint in the bracket positions of this formula was made in || but the right formula was used for calculations. 

The obviously important well known factor exp(2i<7i) written in the equation below was not included in Q. 

Here, the misprint made in Q] is corrected: the factor kq is added in the denominator in front of the square bracket. 

The accurate E Tt , B results for the variants No 3, 4 are given in [li| f in fisll the values of fc res were given by mistake instead 

of -Eres for these variants). Besides, for the case of variant No 14 in [l9t ] the values of a p dC2, a p aC4 are wrong, they should 

be replaced by the values given for this variant in the table of the paper pl| (the lower line). 

There is a misprint in the analogous formula (38) in [j|: the exponent of power in (Nw)~ 2 is absent. 

For the Hulthen potential in Table [I] of the present work we give more accurate calculation results of the ANC and VC for 
3 H and 3 He. 



